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ABSTRACT 

We present a spectral and temporal analysis of the Cathedral QPO detected in the power density 
spectra of the black hole binary and microquasar XTE J1859+226 obtained with RXTE. This peculiar 
type of QPO has been seen on two occasions (on MJDs 51574.43 and 51575.43) during the 1998 outburst 
of this source. It manifests as two peaks with similar amplitudes (~ 3% and ^ 5% RMS) and harmonically 
related centroid frequencies 3 and 6 Hz). The temporal properties of these two peaks are different: 
the amplitude of the ~ 3 Hz feature varies, in anticorrelation with the count rate, by about 50%. The 
~ 6 Hz feature, on the other hand, shows a slight increase 7%) of its amplitude with count rate. 
The RMS-spectra of the two peaks are also quite different. The ~ 3 Hz feature is softer than the other 
one, and although its RMS amplitude increases with energy it shows a cut-off at an energy of ~ 6 keV. 
The RMS of the 6 Hz increases up to at least 20 keV. We also study the bicoherence, i^) of both 
observations. At the diagonal position of the peaks the values &^(~ 3, ^ 3) and 6, ~ 6) are rather 
high and similar to what has been reported in the case of the type C QPOs of GRS 1915-1-105. By 
comparison with the latter source the fact that the bicoherence of the 3 Hz feature is higher than than 
of the other peak, would tend to indicate that the '-^ 3 Hz is the fundamental QPO, and the other its first 
harmonic. The value of 3, ^ 6) is, however, low and therefore indicates a behavior that is different 
than that seen in GRS 1915+105. Wc discuss the implications of these differences in the context of an 
harmonic relationship of the peaks, and suggest that, rather than pure harmonics, we may see different 
modes of the same underlying phenomenon competing to produce QPOs at different frequencies. 

Subject headings: accretion, accretion disks — black hole physics 
J1859-I-226, XTE J1550-564, GRS 1915-1-105) - 
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1. INTRODUCTION 



I Wagoner et al.l 12001, and ref. th erein), Lense-Thirring 



Black Hole Binaries (BHB, e.g.|Remillard fc McClintojg P-eession (e-^^^^ therein) os- 



dilating shocks (e.g. Molteni et al. 1996. and ref. therein), 
or MHD instability (jTaggcr fc Pe llat 1999}). From the ob- 
servational side it seems clear that the inner disk some- 
how sets the fre quency o f LFQPOs (e.g. Muno et al.lll999l : 
[Rodriguez eral . 2002a.bl. l2004bHMikles et al.ll2009D LFO- 
POs are, however, strong when a strong hard component 
is seen in the energy spectra and their fre quency is corre- 



[2001 for a review) transit through different 'spectral 
states' during their outbursts. These are defined both 
by the spectral and temporal parameters respectively ob- 
tained through analyses of their energy spectra and power 
density spectra (PDS). In the soft state, where the emis- 
sion is dominated by a bright and warm (-^ 1 keV) ac- 
cretion disk, the level of variability is weak, and the PDS 
is power law like. On the other hand, in the hard state, 
when the disk is much colder (< 0.5 keV) and thought 
to be truncated at a large distance from the accretor, 
the level of variability is much higher and shows a band 
limited noise component. Other states exist and can be 
seen as intermediate between the Hard and Soft ones. See 
e.g. iRemillard fc McClintockl (|2006l) : iHoman fc Bellonil 
(j2005[ ) for recent reviews and precise classifications of the 
BHB's states. Furthermore, quasi-periodic oscillations in 
the 'low frequency' range (0.1-20 Hz, hereafter referred to 
as LFQPO) are seen in the hard and during the interme- 
diate states. These LFQPOs have been further classified 
into types A, B, or C based on their typ ical frequencies, 
total RMS amplitude, an d time lags (e.g. IRemillard et al.l 
120021 : ICasella et al.l |2005() . We have recently proposed a 
tentative classification of sta tes based on the pre sence of 
the different types of QPOs (|Varniere et al.ll20lH) . 

The exact origin of QPOs is still a matter of debate: 
they could be due to global oscillations of the disk (e.g. 
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France. 



lated with the power law photon index (e .g. lVignarca et al.l 
l2003t IShaposhnikov fc Titarcluig |2007^ . This last point 
would tend to indicate an origin in relation to the corona. 
The amplitude of QPOs increase s with the energy up to 
a cut -off whose energy is variable ([Rodriguez et al.ll2004al 
I2OOI . 

These LFQPOs usually manifest in the PDSs as pow- 
erful peak referred to as the fundamental and a series 
of (sub) harmonics usually of much fainter amplitudes. 
Very little focus has been put so far on these harmonics, 
even if it is clear that they do not completely share the 
same properties as the fundamental. They can show dif- 
ferent s igns for their time lags, as in the case of type B 
QPOs (jCasella et al.112004 ) and show different shapes of 
theRMS spectra (lCui|[l99l iHoman et aLll200lt iRao et al.l 
|2010() . The true identificati on of the fundam ental has re- 
cently been questioned (e.g. lRao et al.|[2"010( ). which raises 
questions related to the genuineness of the harmonic rela- 
tionship of the peaks. 
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XTE J1859+226 was disco vered on 1999 Octo ber 9 with 
the RXTE All Sky Monitor (jWood et al.lll999t) as it was 
entering into outburst. It is a microquasa r given the obser- 
vatio ns of relativistic ejections in radio (jBrocksopp et al.l 
l2002t ). an d the obse r vation s of low and high frequency 
QPOs led ICui et alj (|2000f ) to classify it candidate 
BHB. An extensive timi n g ana lysis of this source is pre- 
sented by ICaseUa et all (|2004[ ). Like XTE J1550-564, 
XTE J1859-h226 displays ah three types of LFQPOs. In 
two particula r obse rvations, during the 1999 outburst, 
ICasella et alJ (|2004[ ) observe the presence of two peaks 
with harmonically related frequencies, but unlike any other 
cases, similar RMS-a mplitudes. These ar e referred to as 
the 'Cathedral' QPO (|Casella et al.ll2004[) : the strongest 
peak (and highest in frequency) has hard lags (the hard 
X-ray lag behind the soft X-rays) and is interpreted as the 
fundamental peak, while the lowest frequency peak has 
soft lags and is the sub-harmonic. Interestingly these two 
observations, although separated by about a day, are also 
separated by an observation showi ng another type of Q PO 
indicative of a spectral transition (jCasella et al.ll2003) . 

The existence of this Cathedral QPOs showing two har- 
monically related peaks of similar amplitude raises chal- 
lenging questions for all theoretical models. To obtain a 
clearer view of the properties of these peaks we performed 
a complete study of the QPO structure from these two ob- 
servations, and compare the temporal and spectral behav- 
iors of the two peaks. The organization of the letter is as 
follows: we start presenting the observations ID and data 
reduction methods. We then describe our results in Sec- 
tion 3, starting with the temporal evolution of the source, 
the fit to the broad band PDS, the energy dependences of 
the QPOs, and we end this part by presenting the bico- 
herence of the observations. We discuss our findings in the 
last section of this letter. 

2. OBSERVATIONS AND DATA REDUCTION 

We focus on RXTE observations 40124-01-24-00 
(Obs. 1) and 40124-01-27-00 (Obs. 2) respectively made on 
1999 October, 23 (MJD 51474.43) and 24 (MJD 51475.43), 
near the peak of the 1999 outburst. The full process of 
RXTE/PCA data reduction was made with the HEASOFT 
v6 . 9 software package. We reduced the Binned and 
Event mode data foll owing standard procedures (see e.g. 
[Rodriguez et al.ll2008[) to obtain light curves filtered from 
low elevation above the Earth, large offset from the source, 
and PCU breakdown. We extracted 7.8125 x 10~^ s binned 
light curves in several energy ranges used for the fine tim- 
ing analysis , and a ^ 2-15 keV light curve from the 
Binned data with a time bin of 8 s (Fig. [T|) to character- 
ize the overall behavior of the source over the observations. 
Note that this range contains most of the counts emitted 
by the source. Power density spectra (PDS) were then 
produced with Powspec vl.O on intervals of 16 s in the 
range 0.0625-64 Hz, all intervals being further averaged 
together before the fitting process. A dynamical PDS 
(DPDS) was also computed between 0.25 and 64 Hz, to 
study the variations of the QPOs with time (Fig. [IJ. 
The PDSs were fitted between 0.0625 and 40 Hz in XSPEC 
V 12. 6.0. The background rate was taken into account 
when estimating the RMS amplitudes of the different 
features following Anet = A^aw x ^5^1 with A the am- 



plitu de, S is the source net rate, and B the background 
rate (iBerger k, van der Klisl[l99l [Rodriguez et al.ll2004al 
120081 ). 

3. RESULTS 

3.1. Time evolution of the source: light curve and 
dynamical PDS 

The 2-15 keV XTE J1859+226 PCA light curve binned 
at 8 s and DPDS of Obs. 1 are reported in Fig. [TJ Large 
variations around a mean raw (net) count rate of 5700 
cts/s (5676 cts/s) are clearly visible. We especially note 
the presence of two rather broad dips lasting respectively 
^ 100 s and ^ 50 s near relative times 1350 s and 1775 s, 
and a third occurring near the end of the observation (at 
t^ 2050 s). The DPDS shows the presence of two strong 
features (relatively to the overall noise) around 3 and 6 Hz. 
These two features are quite thin and indicate the pres- 
ence of QPOs at these frequencies. Interestingly the two 
QPOs seem to have different behavior with time (Fig. [T]): 
the QPO with the smallest frequency is, on average, much 
weaker than the other one, and is strong only when the 
count rate is around its mean value. It is, in particu- 
lar quite weak during the small flares, and is not visible 
during the three dips. The highest frequency QPO, on 
the other hand, seems, in term of power, more stable and 
seems to vary significantly only during the dips (Fig. [T]). 
The feature seems rather broad around 6 Hz, which 
may indicate some rapid variations of the frequency, or 
simply an intrinsically low coherence QPO. A very similar 
behavior (not shown) is also seen in both the light curve 
and DPDS of Obs. 2. 



3.2. Broad band PDS 

We started with fitting the large band PDSs. We did not 
subtract the white noise and preferred adding a constant 
to our fit model t o acc ount for this component. Follow- 
ing ICasella et all (120041) we fitted the continuum with the 
sum of three broad and three narrow Lorcntzians (Fig. [2|) 
on top of the white noise component. This model yields 
a good fit with x^=l-13 (resp. 1.09) for 118 degrees of 
freedom (DOF) for Obs. 1 (resp. Obs. 2). The param- 
eters of the three thin peaks are reported in Table [1] for 
both observations. The parameters of all three peaks are 
compatible between the two observations, which lends cre- 
dence to their complete similarity. In the remainder of the 
paper only the results of Obs. 1 are precisely described 
and presented in the figures. Note that in all cases the 
same analysis was performed on Obs. 2 and the results 
and trends observed are consistent with those of Obs. 1. 
The third peak is compatible with 4 x 1^1, and 2 x j/2 at 
the ~ 2cr level, and is likely to be an harmonic of one of 
the two main peaks. In the remainder of this paper, we 
will refer to either peak 1, 2, or 3, or QPO 1, 2, or 3 for 
the peaks at ~ 2.9, ~ 5.8, and ~ 11.2 Hz respectively. 

Although the fit statistic is quite good, we remark some 
residuals on the lo wer shoulder o f QPO 2. This effect is 
also mentioned by ICasella et al.l (I2004D in Obs. 2, and 
these authors added a Gaussian to better represent the 
peak. In our case, the addition of another thin Lorentzian 
at ~ 5.5 Hz improves the fit to x^=0.90 for 137 DOF, and 
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corrects the defect previously present in the residuals. The 
feature is, however, poorly defined, and its parameters are 
badly constrained. We verify, by re-doing the whole anal- 
ysis that it had no significant impact on the other peaks, 
and since no influence was found it was omitted from our 
study, and is not further discussed here. 

As mentioned in the previous section, QPOi seems more 
intermittent that the second peak (Fig. [Ij. In order to 
quantify and study any possible dependence of the QPOs 
amplitudes with the count rate (Sec. 13. ip. we fol l owed a 
procedure similar to that presented in Heil et al] ()2011[ ). 
Each observation was separated into ten count rate inter- 
vals of equal width. Each interval therefore has a width 
equal to ^^°^('^-^)^~^-^'"('^-^) cts/s. Due to short accumula- 
tion times, and thus poor statistical quality of the resulting 
PDSs, some of pairs of intervals were combined together. 
In the case of Obs. 1 this rebin concerned intervals 1&2, 
and intervals 3&4. The final division of Obs. 1 is repre- 
sented in Fig. [TJ We then extracted a 7.8125 x 10"^ s 
binned light curve per time interval of 8 s, and produced 
a PDS from each of these light curves. Individual PDSs 
belonging to the same count rate interval were averaged 
together to produce the final count rate dependent PDSs. 
The latter were then fitted to estimate the parameters 
of the peaks, with a special focus on their rms ampli- 
tudes. No specific trend is seen between the frequency 
or the coherence of any of the two peaks with the count 
rate. In order to ease the compa rison with the QPOs of 
XTE J1550-564 (jHeil et al.l20lH) . we represented the evo- 
lution of the absolute RMSs (in terms of cts/s) of the two 
peaks with the count rate in Fig.[3]for the particular case of 
Obs. 1. The amplitude of QPOi decreases with increasing 
count rate, while that of QPO2 increases slightly (Fig. [S]) 
and show a possible linear trend. Between the first real 
detections of the two peaks (in the second count rate bin) 
and the last bin, QPOi varies from an RMS amplitude of 
2.7 ± 0.3% to 1.35 ± 0.34 (a variation of 50% in ampli- 
tude) while QPO2 varies from 3.6 ± 0.4% to 3.9 ± 0.2% (a 
non-significant variation of ^ 7.7% in amplitude). 

3.3. Energy dependence of the QPOs 

To study the energetic dependence of the QPOs and 
produce the RMS-Spectra (Fig. d]), we fitted each energy 
dependent PDS with the statistically required number of 
broad features to account for the continuum. We remark 
here that these features also have complex energetic de- 
pendences, and thus the diflFerent fits either require 1,2, or 
3 broad Lorentzians. The study of these is. however, be- 
yond the scope of this letter and we will not discuss them 
further. We then added the thin Lorentzians to account 
for the QPO. The third peak is most of time undetectable 
and no spectrum can be acquired for it. The resultant 
RMS-Spectra for QPOi and QPO2 from Obs. 1 are re- 
ported in Fig. 31 An almost identical behavior is seen in 
the QPO spectra obtained from Obs. 2 (not shown). 

The two QPOs share a common trend: their amplitude 
first increases with energy before reaching a plateau. Such 
QPO-spectra are co mmon for all types (A, B, C) of LFQ- 
POs in this source (jCasella et al.ll200^ . and have, for ex- 
ample , also b een seen in GRS 1915-H 05 (iRodrieuez et all 
2004al [2OOI, and XT E J1550-564 (iHoman et al.ll2001t 
Rodriguez et al.ll2004b[ ). The fact that the normalization 



of the spectra is different is not unexpected since the two 
peaks have different total amplitudes. The precise shape 
and typical parameters (energy of the break, slope, ...) 
are, however, clearly different (Fig. U). This is illustrated 
in the right panel of Fig. 21 where we normalized the RMS 
spectrum of QPO2 by that of QPOi. This representation 
also clearly shows that QPO2 has a steeper (harder) spec- 
trum than QPOi. The latter first increases up to ^ 5.7 
kcV and is then flat until ~ 20 keV. QPO2 is undetectable 
in the first energy bin (and is thus fainter than QPOi). It 
increases up to ~ 20 keV where its plateau is reached. 

3.3.1. Bicoherence of the Cathedral QPO 

The bispcctrum and the related bicoherence permit any 
possible coupling between the different components of the 
PDS to be studied, and thus allow one to go beyond the 
diagn osis brought by the PDSs (see e.g. iMaccarone et al.l 
I2OIII and references therein for a deeper discussion on 
these aspects). Dividing the light curve in K segments of 
equal length, the bicoherence is expressed as: 



where Xi{i') is the v component of the discrete Fourier 



transform of the i*'^ interval (see e.g. Maccarone &: Coppil 



120021 : lUttlevet al.ll2005l : IMaccarone et al.ll201ll for the ide 

scription, restriction, interpretation and applications to 
various astrophysical sources, including BHBs). We cal- 
culated the 2-15 keV bicoherences of XTE J1859-F226 be- 
tween 6.25 X 10~^ Hz and 16 Hz, with a frequency resolu- 
tion of 0.0625 Hz. Fig.jSlshows a zoom on the bicoherence 
plot we obtained in the case of Obs. 1 (the plots and re- 
sults obtained in the case of Obs. 2 are almost identical). 
Note that by definition the bicoherence is symmetric with 
respect to the first diagonal {h'^{v,ii) = b'^{^,v)). The 
highest values of the bicoherence are reached in a region 
around (^1 , vi) (the latter being indicated by the white box 
close to the bottom left corner in Fig. between 2.8 
and ^ 3.2 Hz along the first diagonal (Fig.|S|), with a slight 
broadening at low frequencies. A local maximum is also 
reached around the frequency of QPO 2 (illustrated by the 
second white box in Fig. [5]). Note that the mean value 
of calculated in the white regions represented in Fig. 1 
is higher around QPOi than around QP02. The mean 
value of b^ over the 0.0625-16 Hz range is 0.008 with a 
an RMS of about 0.008. The value of 5^ around (2^2, ^^i) 
(black box in Fig. [5]) is quite low. At the position of the 
two peaks we obtain 6^(1^2, ^^i) ~ 0.01, a value within the 
statistical fluctuations around the mean of 6^, and clearly 
compatible with "noise". Note that the same results are 
obtained when looking at the mean value of b^ in broad 
(e.g. 5x5 frequency-pixels) regions centered on (1^2, 1'l)- 
Overall, this tends to show that either the two peaks are 
not coupled or that their coupling is very weak. 

4. DISCUSSION AND CONCLUSIONS 

We presented here an analysis of some of the prop erties 
of the pecuhar 'Cathedral' QPO (ICaseha et al.ll2004[) seen 
in XTE J1859-I-226. Although we mainly focused here 
on Obs. 1 to illustrate our analysis, similar results and 
trends were obtained in the case of Obs. 2, allowing us 
to use these two observations in our argumentation. The 
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cathedral QPO manifests as two apparently harmonically- 
related peaks of sim ilar amplitudes in the PDSs (Fig. [2]). 
ICasella et al.l (|2004D concluded that QPO2 is the funda- 
mental feature while QPOi is its 'sub'-harmonic. With 
this conclusion, and adding the fact that QPO2 shows 
hard time lags while QPOi show soft time lags, they clas - 
sify the Cathedrals as type B QPOs (jCasella et al.ll2004[) . 



Although QPOi and QPO2 seems to be harmonically re- 
lated, and have been classified as harmonics, their overall 
behavior is quite different. Wc summarize here the main 
differences: 

1. Their amplitude has a different temporal evolution. 
QPO2 is present during the whole observation 
(apart from the 3 dips mentioned in Sec. 13. ip . 
while QPOi is only intermittently seen. Our 
analysis shows that it undergoes larger variations 
of its amplitude than QPO2 (Sec. 13. ip . 

2. The RMS-spectra of these two peaks arc clearly 
different in shape, normalization and typical pa- 
rameters (Fig. In particular, the characteristic 
energy at which their spectra flatten is different by 
about a factor of > 3: QPOi reaches its peak at 
around ~ 6 keV, while QPO2 is much harder and 
peaks at > 20 keV. 

3. The time lags of the two peaks are different: QPOi 
has soft lags (the soft X-rays lag behind the ha rd 
ones) , while QPO2 has hard lag (jCasella et al.l 
120041) . Thi s property is a definition of type-B 
QPO s (e.g. iRemillard et"all 120021: ICasella et all 
[2001 . 

4. The value of the bicoherence h^{v2, 1^1) is compati- 
ble with statistical noise, which tends to indicate a 
very weak or no coupling between the peaks at 1/2, 
vi , and vi + V2. 

While poi nts 2, and 3 have been mentioned for other 
BHBs fc.g. iHoman et al.ll2001l: iRao et al.l[20Tot [Cuilll999l 

for XTE J1550-5 64 and GRS 1915-| -105), in addition to 
XTE J1859+226 (jCasella et al.ll2004[) . it is the first time 
to our knowledge that points 1 and 4 are rep orted for 
any BHBs. In XTE J1550-564 rHeil et al.l ()2011h report a 
positive linear RMS-flux relation for (type-C) QPOs with 
a frequency smaller than ~ 4 Hz, and a negative RMS-flux 
trend at higher frequencies. Here the situation is opposite. 
QPO2, that has the highest frequency, shows the linear re- 
lation, and QPOi shows the negative trend (Fig. [3]). This 
could be the signature of a new fundamental difference 
between types B and C QPOs. 

In XTE J1550-564, iRao et al.l (|2010[ ) mentioned that 
the difference of energy spectra of the two peaks, especially 
the fact that the fundamental has a higher amplitude at 
high energy can indicate a more sinusoidal signal at higher 
energies. In this respect, the fact that QPOi is stronger at 
low flux could also indicate a more sinusoidal signal in the 
peaks. This explanation, although simple and tempting, 
do not account for the different signs of the time lags of the 
two peaks, since one would naively expect in the case of a 
non-purely sinusoidal signal that all components undergo 
the same physical processes and thus show similar lags. 



In addition, iRao et aTl (j2010t) come to the conclusion that 
the peak referred to as sub-harmonic in XTE J1550— 564 
could in fact be the true fundamental feature. In this 
case the "more sinusoidal" explanation does not hold any 
more since the harmonic (and therefore the peak with the 
highest frequency) should then disappear first, and show 
a softer spectra, which is clearly opposite to what is seen 
here. 

It is interesting to remark, here, that the bicoherence 
behavior of the two peaks (at {vi^vi) and {1^2,1^2)) is 
s imilar to that of the type C QPOs of GRS 1915-^105 
CMaccarone et all |2011[ ) , where reaches a high value 
close to the frequencies of the peaks. In this latter 
source is high at the (bi-)position of the fundamen- 
tal, a nd much lower for the harmonic (jMaccarone et all 
l2011t ). Pursuing the comparison with GRS 1915-1-105, in 
XTE J1859-H226 the values of fo^ would tend to indicate 
that QPOi is the fundamental and QPO2 its harmonic. 
In this latter source, however, different global patterns 
have been identified in the bicohrence plot. In all cases 
(where a strong harmonic is seen) a strong coupling is 
seen between the fundament al and the harmonic (e. g. the 
"web" or "cross" patterns accarone et all 120 111) , and 
even between the noise component and the QPO. This 
is not observe in the case of XTE J1859-I-226 (Fig. rcf- 
fig:bicoh). 

An easy interpretation for the four points above would 
be that the two peaks have absolutely no relation together 
and that the apparent harmonic relation of their frequen- 
cies is fortuitous. This seems difficult to reconcile with the 
fact that a similar behavior is seen both in Obs. 1 and 2 
that are separ ated by roughly a d ay, and a transition into 
another state (ICasella et al.|[200^ . In addition, point 2 in 
the case of type C QPOs, and points 2 and 3 in the case of 
type B QPOs have been reported i n this and other sources 
Cuil I1999I: iHoman et al l I2OOII: iRemiUard etlH 120021: 



Rodriguez et al.l l2002al : ICasella et al.l 12004 iRao etal l 
20101) . It is more likely that the two peaks are somehow 



related or that, at least, a possible common mechanism 
sets their frequencies to be integer multiple, with the two 
features not being harmonics in a physical sense. 

Hard lags are usually easy to understand in the context 
of Comptonization of soft photons in a h ot and ten uous 
medium (e. g. the so-c alled corona; e.g. |Cui"l99?). As 
discussed inlCuH (|1999D . and lGierlihski fc Zd^iarski (2(305l) 
the energy dependences of both the QPO and/or the con- 
tinuum are indicative of variations of physical parameters: 
the favored ones are eith er the tem perature and/or op- 
tical depth of the corona (|Cuilll999l ). the temperature of 
the soft pho tons, or a variable power of the Comptonized 
component (jGierlihski fc Zdziarskill2005l ). In this respect, 
however, the different time dependences of the two peaks, 
and the opposite sign of their time lags are difficult to 
understand. 

Other families of models involve precession at the inner 
boundary of the disk, but here again, the fundamental and 
harmonics are produced at the same location and should 
'sec' the same environment. In the framework of Comp- 
tonization the time lags, RMS-spectra, and bicoherence of 
the different features should be the same. 

A way to possibly reconcile the four points summarized 
above could be that the peaks are not harmonics in the 
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usual sense but that they represent different modes of the 
same mechanism favored at different moments, depending 
of some (unknown) parameter (s) in the corona-disk-jet 
system. This hypothesis has the advantage of providing 
an explanation for the time behavior of the two peaks, and 
permits in particular that they do not necessarily appear 
at the same time, but enter in a sort of competition. The 
competition has also been m entioned for other t ypes of 
QPOs in XTE J1859-f226 bv lCasella et al.l (I2004D and in 
the case of XTE J1550-564 by iRemillard et al.l (|2002| ) 
although it has never been explored further. The physical 
states of the disk/corona/jct system might then set the 
conditions for one or the other peak to dominate. This will 
be explored by us from a particular theoretical standpoint 
in a forthcoming paper (Varniere, Tagger, Rodriguez, sub- 
mitted to A&A). we note that the behavior observed in 
the case of the cathedral QPO is different from, at least, 
that of the type C QPO of GRS 1915-^105 with respect 
to the RMS-flux relation and to the bico herence behavior 
(jHeil et al.ll201ll: iMaccarone et al.ll201lt) . This may, how- 
ever, represent a fundamental property of type B QPOs 



in general. In any cases, the new mo del- independent ob- 
servational facts presented here should be taken as strong 
constraints in any attempts to model LFQPOs in BHBs. 
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Table 1 

QPO PARAMETERS OBTAINED FROM THE FITS TO THE 2-40 KEV PDS. THE ERRORS ARE GIVEN AT THE 90% LEVEL. 

Obs. # VI Qi* Ai^ V2 Q2* Qs* 

(Hz) (% RMS) (Hz) (% RMS) (Hz) (% RMS) 

1 2.94 ±0.02 5^9 2.8 ±0.1 5.828 ± 0.025 7^3 4.7 ±0.1 11. 2l^-^ 9^5 1.1+^-^ 

2 3.00 ±0.04 5.2 2.9 ± 0.2 5.86 ± 0.04 6.5 4.6+g:i 11.5 ±0.3 8.3 1.5to:3 



* Q=!//FWHM. 

A stands for RMS amplitude. 
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Fig. 1. — Upper panel: 2-40 kcV dynamical power spectrum of XTE J1859+226. Lower panel: 2-40 keV light curve. 
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Fig. 2. — 2-40 koV power spectra showing the presence of the 'cathedral' QPO, the best fit model, and the single components used in the 
fit. The solid lines show the thin features (QPO and harmonic), the dashed line show the continuum component (including the white noise). 




Frequency (Hz) 



Fig. 3. — Evolution of the QPOs absolute amplitudes with the count rate for Obs. 1. The diamonds correspond to QPOi and the circles 
to QPO2. 
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The cathedral QPO in XTE J1859+226 
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Fig. 4. — Left: Energy dependences of the two main QPOs. Upper limits are given at the 90% level. Right: Ratio between the 
amplitudes of QPO2 and QPOi. 
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Fig. 5. — Bicoherence plot for Obs. 1 on a frequency region including the two peaks. The white regions are 5x5 pixels boxes centered on 
the positions of the two peaks (at respectively ~ (2.94, 2.94) Hz and ~ (5.83, 5.83) Hz), while the black one is centered on ~ (5.83, 2.94) Hz 
and shows the value of at the intersection of the two peaks. 

Frequency (Hz) 




0.05 0.068 0.086 0.1 0.12 



